Localization of a Bose-Einstein condensate in a bichromatic optical lattice 
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By direct numerical simulation of the time-dependent Gross-Pitaevskii equation we study different 
aspects of the localization of a non-interacting ideal Bose-Einstein condensate (BEC) in a one- 
dimensional bichromatic quasi-periodic optical-lattice potential. Such a quasi-periodic potential, 
used in a recent experiment on the localization of a BEC [Roati et al, Nature 453, 895 (2008)], can be 
formed by the superposition of two standing-wave polarized laser beams with different wavelengths. 
We investigate the effect of the variation of optical amplitudes and wavelengths on the localization 
of a non-interacting BEC. We also simulate the non-linear dynamics when a harmonically trapped 
BEC is suddenly released into a quasi-periodic potential, as done experimentally in a laser speckle 
potential [Billy et al, Nature 453, 891 (2008)]. We finally study the destruction of the localization 
in an interacting BEC due to the repulsion generated by a positive scattering length between the 
bosonic atoms. 
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I. INTRODUCTION 



The localization of the electronic wave function in 
a disordered potential was predicted by Anderson fifty 
years ago [![. More recently the phenomenon of local- 
ization due to disorder was experimentally observed in 
electromagnetic waves [U, H| , in sound waves [1] , and also 
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in quantum matter waves 0, @, 
quantum matter waves, Billy et al. 
tial localization of a Bose-Einstein condensate (BEC) of 
87 Rb atoms released into a one-dimensional (ID) waveg- 
uide in the presence of a controlled disorder created by a 
laser speckle. Roati et al. 0] observed localization of a 
non-interacting BEC of 39 K atoms in a ID potential cre- 
ated by two optical-lattice (OL) potentials with different 
amplitudes and wavelengths. The non-interacting BEC 
of 39 K atoms was created by tuning the inter-atomic 
scattering length to zero near a Feshbach resonance Q. 

In this paper, with intensive numerical simulations of 
the linear Gross-Pitaevskii (CP) equation, (which is just 
the Schrodinger equation,) we study different aspects 
of localization of the BEC in a ID bichromatic quasi- 
periodic OL potential used in the experiment of Roati et 
al. [f|. The ID quasi-periodic potentials have a spatial 
ordering that is intermediate between periodicity and dis- 
order [l(J HH, G3 • In particular, the ID discrete Aubry- 
Andre model of quasi-periodic confinement [ll|, G3 dis- 
plays a transition from extended to localized states which 
resembles the Anderson localization of random systems 
[HI, [HI • Modugno [HI has recently shown that the lin- 
ear ID Schrodinger equation with a bichromatic periodic 
potential can be mapped in the Aubry- Andre model and 
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he studied the transition to localization as a function of 
the parameters of the periodic potential. 

To investigate the interplay between the bichromatic 
potential and the inter-atomic interaction in the localiza- 
tion of a BEC, we adopt the ID non-linear GP equation 
[r| [13] in place of the linear ID Schrodinger equation 
used to describe a non- interacting BEC. We find that 
the non-linearity of the GP equation, which accurately 
models the binary inter-atomic interaction of atoms, has 
a strong effect on localization and a reasonably weak re- 
pulsive non-linear term is capable of destroying the lo- 
calization. Our results on the effect of non-linearity in 
the localization are thus in qualitative agreement with 
similar predictions based on the ID discrete non-linear 
Schro ding er equation (DNLSE) with random on-site en- 
ergies |18| . First effects of a weak non-linearity in Ander- 
son localization have been shown experimentally in light 
waves in photonic crystals [1, [lj| . 

There have already been a number of theoretical stud- 
ies on Anderson localization under the action of different 
potentials. Sanchez-Palencia et al. and Clement et al. 
considered Anderson localization in a random potential 
fl3| . Damski et al. and Schulte et al. considered Ander- 
son localization in disordered OL potential [3]. There 
have been studies of Anderson localization with other 
types of disorder [2(j. Effect of interaction on Anderson 
localization was also studied [2l|. Anderson localization 
in BEC under the action of a disordered potential in two 
and three dimensions has also been investigated [13| ■ In 
this paper we study different aspects of the localization 
of an ideal BEC and the delocalization of an interacting 
BEC in a quasi-periodic OL potential using the linear 
and non-linear ID GP equation. 

In Sec. [II] we present a brief account of the non-linear 
ID GP equation used in our study and of the variational 
solution of the same under appropriate conditions. In 
Sec. IIIII we present our numerical studies on localiza- 
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tion using time propagation under the Crank-Nicolson 
discretization scheme. The density profile of the tightly 
localized states are in agreement with the variational re- 
sults. We also study the effect of the variation of the 
wavelength and intensity of the OL potentials on the lo- 
calization. We study the non-equilibrium dynamics, as 
observed in the experiment of Billy et al. [5(, when a 
harmonically trapped BEC is suddenly released from the 
harmonic trap into a quasi-periodic OL potential. We 
also investigate the destruction of localization in the pres- 
ence of a repulsive atomic interaction and it is found that 
a reasonably small non-linearity can destroy the localiza- 
tion of the BEC. In Sec. HVI we present a brief discussion 
and concluding remarks. 



II. ANALYTICAL CONSIDERATION OF 
LOCALIZATION 

In the actual experiment of Roati et al. [6| the ID 
quasi-periodic bichromatic OL potential was produced 
by superposing two OL potentials generated by two 
standing-wave polarized laser beams of slightly different 
wavelengths and amplitudes, which we take here simi- 
lar to those in the experiment of Roati et al., e.g., with 
wavelengths Ai = 1032 nm and A2 = 862 nm. This ID 
quasi-periodic OL potential can be written as Q 



lowing GP equation [16|, [lj 



V(z) 



2siEi cos 2 (fcjz), 



(1) 



where 2s,-,i = 1,2, arc the amplitudes of the OL 
potentials in units of respective recoil energies Ei = 
2tt 2 H 2 /(mAj), and fcj = 27r/A;, i = 1,2 are the respective 
wave numbers, Ti is the reduced Planck constant, and m 
the mass of an atom. 

With a single periodic potential of the form cos 2 (kz) 
with S2 = 0, the linear Schrodingcr equation permits only 
de-localized states in the form of Bloch waves. Localiza- 
tion is possible in the linear Schrodingcr equation due to 
the "disorder" introduced through a second component 
in Eq. (JTJ) . The localization is not intuitively obvious. 
The potential ([1]) continue to have an infinite number of 
finite barriers as in a simple OL potential, and it might 
be expected that any localized state will decay due to 
tunneling. 

The localized states that we study are low-lying states 
of the system with potential (fTJ). They are quite distinct 
from the so called gap solitons in a simple OL poten- 
tial appearing for repulsive non-linearity in the band -gap 
of the spectrum of the linear Schrodingcr equation [231 ] - 
These gap solitons with finite spatial extension are ex- 
cited states of the system without a linear counterpart. 

The dynamics of a trapped BEC of N atoms in a trans- 
verse harmonic potential of frequency plus the axial 
quasi-periodic OL potential |T]) is determined by the fol- 
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with normalization j \ip(r, t)\ 2 dr = 1 and where g = 
2Na/a± is the dimensionless interaction strength with a 
the inter-atomic scattering length and aj_ = \ l /K/(mu)±) 
the characteristic harmonic length of the transverse con- 
finement, and r = (x,y,z) defines the Cartesian coor- 
dinates. Here we have considered the harmonic trap 
(x 2 + y 2 )/2 in transverse directions (x, y) and the quasi- 
periodic potential V{z) in the longitudinal direction z. 
In Eq. (|2|) length is in units of a± , time in units of wj_ , 
and energy in units of Hlu±. In the non- interacting case 
g = and Eq. ([2]) becomes the usual linear Schrodingcr 
equation. 

Another completely equivalent potential is the one 
where the cosine term of Eq. (|3|) is replaced by a sine: 
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However, potential ([4]) generates a different type of local- 
ized states compared to potential ((3]). Potential j4} has 
a local minimum at the center z = 0, consequently stable 
stationary solutions with this potential have a maximum 
at z = 0. However, potential ([3]) has a local maximum 
at z = corresponding to a minimum of the stationary 
solution at the center. We shall show that, starting with 
an initial Gaussian wave function centered at z = 0, the 
numerical solution of the GP equation ((2]) gives differ- 
ent localized cigenstates depending on the choice of the 
confining potential. 

For a cigar-shaped trap with strong transverse con- 
finement, it is appropriate to consider a ID reduction of 
Eq. @ by freezing the transverse dynamics to the re- 
spective ground state and integrating over the transverse 
variables. The resulting 3D-1D reduction of Eq. for 
small non-linearity g is realized via [24|, [25| 
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(z,t), (5) 



with normalization dz\<fi(z, t)\ 2 = 1 and where (f>(z, t) 
is the axial wave function of the Bose condensate. Eq. 
([5]) is a ID non-linear Schrodingcr equation with cubic 
non-linearity. 

Modugno used potential (j4|) in his study and in addi- 
tion took Ai = 2tt and defined = 2sj, and (3 = A1/A2 = 
2tt/A 2 , so that (viz. Eq. (16) of Ref. (H) 
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Eq. ((5|) with potential ((6|) can be mapped flal in a non- 
linear version of the Aubry- Andre model [id . by ex- 
panding the wave function (j>(z, t) = J^j c j (t)Wj (z) over a 
set of orthonormal Wannier states Wj(z), where Wj(z) is 
maximally localized at the j-th minimum of the primary 
lattice. In this way one finds that the complex coeffi- 
cients Cj(t) satisfy the time-dependent discrete non-linear 
Schrodinger equation [l5| 

^ c j(*) = - J ( c i+i(*) + c i-iO)) 

+ Acos(2n/3j)c j (t) + g\c j (t)\ 2 c j (t) (7) 

where \cj (t)\ 2 gives the probability of occupation of the j- 
th site at time t. The hopping term J and the "disorder" 
term A arc connected to the parameters of the bichro- 
matic potential ([1]) and can be calculated by using the 
Wannier functions [3, Hfjj]. Modugno [HI showed that 
the linear (g = 0) GP equation J5]) displays localization 
for a large enough A/ J and the localization increases as 
A/ J increases for a wide range of values of /3. 

The Fourier transformation f(k) of the Anderson lo- 
calized state <f>{z) is defined as 

/oo 
exp^Trfcz)^), (8) 
-00 

and the momentum distribution of the localized state 

p(k)^p(k). 

Usually the stationary bound states formed with quasi- 
periodic OL potentials ([3]) and Q occupy many sites 
of the quasi-periodic OL potential and have many local 
maxima and minima. For certain values of the param- 
eters potential (0| leads to bound states confined prac- 
tically to the central site of the quasi-periodic OL po- 
tential. When this happens, a variational approximation 
with Gaussian ansatz leads to a reasonable prediction for 
the bound state. 

The stationary form of the linear Schrodinger equa- 
tion JS]) (with id/dt replaced by a chemical potential fi) 
with potential ^ can be derived from the following La- 
grangian 

/oo 
H 2 (z) - W(z)} 2 /2 - V{z)4> 2 {z)] dz - /i, (9) 
-00 

by demanding 6L/5(f> = 8L/5)jl = 0, where \i is the prime 
denotes space derivative. To apply the variational ap- 
proximation we use the Gaussian ansatz [27] ] 

^ ) = vr -V4^ exp f__^y (10) 

where variational parameters are the norm A/", width w, 
and \x. This ansatz implies that the center of the sta- 
tionary state is placed at the local minimum at z = of 
the quasi-periodic OL potential (H}. The substitution of 
ansatz (|T0|) in Lagrangian ([9]) leads to 

2 

L = fj,(Af—l) — +J2 M[exp(-a> 2 ) - 1], (11) 



where Ai = Aii 2 Si/\ 2 , a, = 27r/Ai. The first variational 
equation from Eq. (TTTj) , dL/d[i = 0, yields Af = 1, which 
will be used in other variational equations. The second 
variational equation dL/dw = 0, yields 

2 

1 = 2otiAiW A exp{-a 2 w 2 ) 7 (12) 
i=i 

and determines the width w. The last variational equa- 
tion dL/dAf=0, yields 

2 

,i = 1/(4 W 2 ) - A t [eM-alw 2 ) - l]/2, (13) 
t=i 

which defines the chemical potential. (However, some 
caution should be exercised in using these variational 
equations: they predict a false bound state for a one- 
term periodic potential with S2 = 0, cf. Eq. (|12p . and 
the one-term potential is known to support no localized 
bound state.) 

III. NUMERICAL RESULTS 

To perform a systematic numerical study of localiza- 
tion with potentials ([3]) and ([4J we vary Ai and A2 main- 
taining the ratio A2/A1 = 0.86 (roughly the same ratio 
A2/A1 as in the experiment of Roati et al. Q). Wc 
consider a transverse harmonic oscillator length oj_ ps 1 
/im, so that the experimental wavelengths Ai = 1032 nm 
and A2 = 862 nm in dimensionless units become Ai ~ 1 
and A 2 « 0.86. In the first part of the study we also 
set scattering length a — corresponding to ideal non- 
interacting bosons. 

We perform the numerical simulation employing 
mostly real-time propagation with Crank-Nicholson dis- 
cretization scheme [28| with space step 0.025, time step 
0.0005. Imaginary-time propagation routine can deter- 
mine the strongly localized state confined to a single OL 
site or so in an efficient fashion. However, imaginary-time 
propagation routine demonstrated difficulty for weakly- 
confined state extending over a large number of OL sites. 

Because of the oscillating nature of the potential great 
care was needed to obtain a precise localized state. The 
accuracy of the numerical simulation was tested by vary- 
ing the space and time steps as well as the total number 
of space steps. A larger value of the ratio s 2 /si gives 
more binding for the localized state, consequently, the 
localized state has smaller spatial extension. 

To understand the nature of these localized states and 
their behavior under the variation of different parame- 
ters, first we consider the localized states with larger val- 
ues of Ai. Such states with large S2/si(= 1) occupy a 
small number of OL sites and hence their simulation can 
be performed relatively easily. The shape of the local- 
ized state then becomes a quasi-Gaussian for potential 
and we compare our numerical results with the vari- 
ational results in this case. In Fig. [1] (a) we plot the 
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FIG. 1: (Color online) (a) Typical density distribution ^?{z) 
vs. 2 for a non- interacting BEC for potential @ for A2/A1 = 
0.86, si = S2 — 2, and Ai = 15, 10 and 5. The numerical 
results are shown by continuous lines. The variational results 
are shown by a chain of symbols, (b) Typical density distri- 
bution 4> 2 {z) vs. 2 for potential ([3]) for same parameters as in 
(a). Both 4> 2 { z ) arm z are in dimensionless units. 



results of density distribution from numerical and vari- 
ational calculations for potential ^ for different Ai for 
fixed si, S2, and A2/A1. In Fig. [T](a) real- and imaginary- 
time propagation routines produced identical results. In 
Fig. [T](b) we plot the density distribution <fi 2 (z) vs. z for 
potential ([3]) for the same set of parameters as in Fig. Q] 
(a). In this case one has two peaks in the density distri- 
bution. We also calculated the energies of these states. 
The energies for Ai = 15, 10, and 5 for the potential (U]) 
of Fig. [I] (a) are 0.264, 0.599 and 2.374 to be compared 
with the variational results of Eq. (fHij) : 0.266, 0.594, and 
2.396, respectively. This agreement between the numeri- 
cal and variational results of density for potential @ in 
Fig. [T] (a), and of the respective energies, provides as- 
surance about the accuracy of the numerical code used 
in simulation in our investigation. We also calculated 
the (numerical) energies of the density profiles displayed 
in Fig. Q] (b) which are, respectively, 0.270, 0.603, and 
2.412. These energies for potential are distinct from 
those of potential ((3]); but the two sets of energies are 
very close to each other. 
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FIG. 2: (Color online) (a) Potential V(z) vs. z given by Eq. 
© for Ai = 2,A 2 /Ai = 0.86, si = 2, and s 2 /si = 0.2 and 
(b) density distribution 4> 2 (z) vs. z of a non-interacting BEC 
calculated with this potential for si = 2, S2/S1 = 0.2, A2/A1 = 
0.86 and Ai = 2,4,8, and 12. (c) Density distribution (j) 2 {z) 
vs. 2 of a non-interacting BEC for potential (T3| for si = 2, 
82/81 = 0.2, A2/A1 = 0.86 and Ai = 2,4,8,12. All variables 
are expressed in dimensionless units. 



Now we present the results for the solution of Eq. 
([5]) with potential ([3]) for a small s^/si. The potential 
V{z) given by © for Ai = 2, A2/A1 = 0.86, si = 2 and 
S2/S1 = 0.2 is plotted in Fig. [5] (a). This potential is 
quite similar to the potential illustrated in Fig. 1(a) of 
Roati et al. @. The density distribution 4> 2 {z) vs. z cor- 
responding to the localization for this potential is plotted 
in Fig. H (b) for si = 2,s 2 /si = 0.2, A2/A1 = 0.86 and 
Ai = 2,4,8, and 12. It is seen that with a decrease of 
Ai more attraction is created. Consequently, the states 
with Ai = 2 and 4 are more localized in space with lesser 
spatial extension. But with further increase of Ai a single 
site of the quasi-periodic OL potential occupies a large 
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FIG. 3: (Color online) Momentum distribution P(k) vs. k 
of the localized states, shown in Figs. (b) and (c), for (a) 
potential ([3| and (b) potential (Q, for A2/A1 = 0.86, si = 
2, sa/si = 0.2, and Ai = 2, 4, 8, 12. 



region in space. When this happens the spatial size of 
the localized bound states increases with Ai , as the local- 
ized state cannot occupy less than two sites of the quasi- 
periodic OL potential (the density for potential © has 
to be symmetric around z = and must have a minimum 
at z = 0). This happens for Ai = 8 and 12. For Ai = 8 
and 12 the localized bound states occupy 4 sites of the 
quasi-periodic OL potential. Partially similar feature is 
also exhibited by the localized states of the quasi-periodic 
OL potential (0| for larger values of Ai as shown in Fig. [2] 
(c), where we plot the density distribution <fr 2 (z) vs. z for 
s 2 = 2,s 2 /si = 0.2,A 2 /Ai = 0.86 and A x = 2,4,8, and 
12. As expected, the localized state with potential ([3]) 
has a minimum at z = 0, whereas the localized state of 
potential (01 has a maximum at z = 0. We plot in Figs. 
|3] (a) and (b) the momentum distribution P(k) vs. k for 
the localized states — a quantity of experimental inter- 
est [f| — shown in Figs. [3] (a) and (b) for potentials © 
and ([?]), respectively. In general, as expected, the central 
peak of momentum distribution of the localized state is 
sharper for a localized state of larger spatial extension. 

Now we study how these localized states are affected 
by a variation of the ratio S2/S1 when Ai,A2 and si 
are maintained constant. For this purpose we plot in 



FIG. 4: (Color online) Typical density distribution 4> 2 (z) vs. 
z for a non-interacting BEC for (a) potential (JJ) and (b) po- 
tential © for Ai = 10, A 2 = 8.6, si = 2, and s 2 /si = 0.2,0.4 
and 1. The quasi-periodic OL potential V(z) for S2/S1 = 0.2 
is also plotted in (a) and (b) in arbitrary units. All variables 
are in dimensionless units. 



Figs. |4] (a) and (b) the density distribution of the lo- 
calized states for potentials (U]) and ©, respectively, for 
Ai = 10, A2 = 8.6, Si = 2 and for different values of the 
fraction s 2 /si = 0.2,0.4 and 1. In addition we plot the 
quasi-periodic OL potential V(z) for 82/ S\ = 0.2 in arbi- 
trary units, just to compare the position of the maxima 
and minima of the potential with the position of the min- 
ima and maxima of density. (The position of the maxima 
and minima of the potential V(z) does not change visi- 
bly with S2/S1.) It is found in both cases that the states 
with larger values of the fraction S2/ s\ have smaller spa- 
tial extension corresponding to larger attraction. The 
results reported in Figs. Q] (a) and (b) arc in qualitative 
agreement with a conclusion of the study of Modugno 
[15(, that the localization appears and increases as the 
disorder to hopping ratio A/ J increases. The increase of 
A/ J is related to an increase of S2/ s\ for a fixed Ai and 
A2, exactly as illustrated in Figs. |U(a) and (b). However, 
because of the distinct model (DNLSE) used in Ref. [HI , 
in contrast to the numerical solution of the GP equation 
in the present study, a quantitative comparison of the 
results of the two studies is not to the point. 
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FIG. 5: (Color online) Typical density distribution 4> 2 {z) vs. 
z for an interacting BEC for different g = 2Na/a± for (a) po- 
tential (0} and (b) potential © for Ai = 4, A2/A1 = 0.86, si = 
2, S2/S1 = 0.2. The quasi-periodic OL potential V(z) is plot- 
ted in arbitrary units, (c) The rms size vs. non-linearity g of 
the stable condensate in the quasi-periodic potential ((3]) (co- 
sine) and (sine) for Ai = 4, A2/A1 = 0.86, si = 2, s^/si = 
0.2. All quantities plotted are in dimensionless units. 



Next we study the effect of including interaction in a 
BEC of 39 K atoms with scattering length a = 33ao = 
1.75 nm [29| (with a = 05292 nm, the Bohr radius,) by 
solving Eq. ([5]) with potentials ([3]) and (|4]). In present 
dimensionless units this will correspond to a scattering 
length of a = a/a± = 0.00175. The inclusion of the repul- 
sive non-linear potential term in Eq. ([5]) will reduce the 



FIG. 6: (Color online) (a) The rms size of an expanding BEC 
released from a harmonic trap into the quasi-periodic OL po- 
tential Q with Ai = 10, A 2 = 8.6, si = 2,s 2 /si = 0.2 at 
time t = 0. The rms sizes during expansions in a periodic 
OL potential with S2 — and free expansion are also shown, 
(b) The density profiles during above expansion in the quasi- 
periodic OL potential are illustrated by plotting <f> 2 (z,t) vs. 
z and t. 



possibility of the appearance of localized bound states. 
This is illustrated in Figs. [5] (a) and (b) where we plot 
the density distribution for Ai = 4, A2/A1 = 0.86, Si = 
2, S2/S1 = 0.2 for potentials (0]) and respectively, for 
different g = 2Na/a± = (0,2,4,5). The corresponding 
quasi-periodic OL potentials are also plotted in arbitrary 
units. (The advantage of using the variable g, rather 
than the scattering length a and number of atoms N, in 
these plots is that the present plots can easily be used 
to simulate different experimental situations with differ- 
ent traps and distinct bosonic atoms.) For g = 0, for 
both potentials the localized states are confined between 
z = ±10. For g = 2, with increased repulsion, the mat- 
ter density is reduced in the central peaks and new peaks 
appear for larger z values. For g = 4, with further in- 
crease in repulsion the matter density is further reduced 
in the central region and new peaks appear in the form 
of ondulating tails near the edges. With further increase 
in the value of g, the localized states have larger and 
larger spatial extension and soon the non-linear repul- 
sion is so large that no localized states are possible and 
this happens rapidly as g is increased beyond 6. The 
non-linearity in Eq. J5]) is g — 2aN/aj_ and for about 
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180 39 K atoms with a = 0.00175 ^ the non-linearity 
has the typical numerical value g ss 6. Such a small 
non-linearity can have a large effect on localization of a 
39 K BEC and can prohibit the localization. However, the 
number of K atoms can be proportionately increased if 
the scattering length is reduced by varying an external 
background magnetic field near a Feshbach resonance Q . 
As g value is increased, the root mean square (rms) size 
of the BEC increases before reaching a critical g value 
(> 6) corresponding to the destruction of localization. 
(It is difficult to obtain accurately the critical value of g 
needed to destroy the localization.) The increase in the 
rms size of the localized state with the increase in g is 
illustrated in Fig. [5] (c) where we plot the rms size vs. 
g for potentials j3]) and ((4]). It should be noted that in 
the experiment of Roati et al. @ the residual scattering 
length of 39 K atoms near the Feshbach resonance was 
O.lao (= 0.0053 nm), e.g., they can vary the scattering 
length in such small steps. Thus it should be possible 
experimentally to obtain the curves illustrated in Fig. [5] 
(c) and compare them with the present investigation. 

So far we studied the stationary properties of the lo- 
calized state. Next we study some dynamical aspects of 
the localization. One interesting feature is what hap- 
pens when a BEC is released from a harmonic trap into 
a quasi-periodic OL trap as investigated in the experi- 
ment of Billy et al. Q . Mathematically, it means that at 
time t = the harmonic trap is suddenly changed into a 
quasi-periodic OL trap. The evolution of the rms size of 
the BEC with time is plotted in Fig. [6] (a) for potential 
(TJ} . The rms size first increases with the expansion of the 
BEC and after a certain amount of expansion it will be 
locked in an appropriate localized state. After this hap- 
pens the system executes breathing oscillation around a 
mean shape of the localized state and the rms size re- 
mains bounded and does not increase indefinitely with 
time. To compare with this behavior we also plotted the 
rms size for expansion in a pure periodic OL potential 
when the disorder is removed by setting s% =0 in poten- 
tial Q. In that case there is no localized state and the 
system expands for ever with an ever-increasing rms size. 
Finally, when all OL potentials are removed by setting 
S\ = S2 = the system increases monotonically with a 
higher expansion rate. The dynamical density profile of 
the BEC during the expansion and locking in a local- 
ized state for the quasi-periodic potential ^ is shown in 
Fig. [6] (b) where the initial expansion until t = 80 and 
the consequent breathing oscillation of the BEC is clearly 
illustrated. 



IV. CONCLUSION 

In this paper, using the numerical solution of the GP 
equation, we studied the localization of a non-interacting 
BEC in a quasi-periodic ID OL potential prepared by two 
overlapping polarized standing-wave laser beams with 
different wavelengths and amplitudes. Specifically, we 
considered two analytical forms (sine and cosine) of the 
OL potential. We studied the effect of the variation of 
wavelengths and amplitudes on the localization. We also 
studied the non-linear dynamics when a BEC is released 
from a ID harmonic trap into a quasi-periodic ID OL 
trap. After release, the BEC first expands (from the 
tightly bound harmonic oscillator bound state) and then 
the expansion is stopped and the BEC is found to be 
trapped into one of the localized states in the quasi- 
periodic OL potential, as observed by Billy et al. [5j. 

We also studied the effect of a repulsive atomic in- 
teraction on localization. It is found that a repulsive 
atomic interaction destroys the localization. We studied 
the route to this destruction of localization in some de- 
tails. In particular, we investigated the localization as 
the non-linearity g = 2Na/a± of the non-linear ID GP 
equation is increased. It is found that as g is gradu- 
ally increased, the localization is slowly weakened with 
the localized state extending over a large space domain. 
Eventually, for g greater than about 6 or so, the localiza- 
tion is destroyed. 

There have been previous studies of some aspects of 
the localization of a BEC in a quasi-periodic potential 
[l5[ and also its destruction [l8| . (It should be noted that 
the present study is mostly complimentary to these pre- 
vious studies, rather than overlapping.) However, in the 
present study we consider a direct numerical solution of 
the GP equation as opposed to a solution of the DNLSE 
used in the previous studies. In view of the rapidly os- 
cillating nature of the quasi-periodic OL potential and of 
the solution of the non-linear equation, the results of the 
direct numerical solution of the GP equation, as used in 
the present investigation, should be more useful for a 
direct comparison with the experiments. 
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